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Natural  electromagnetic  fields  at  the  floor  of  the  deep  ocean  have  been 
us«d  to  infer  the  electrical  conductivity  of  the  upper  mantle  using  the  pas¬ 
sive  magnetotellur ic  method  (Coz  sA  al .  1970,  Fillouz  1980,  Chave  e£  ml. 
1981).  All  such  studies  indicate  rising  electrical  conductivity  at  depths  of 
50-100  km,  but  the  data  cannot  resolve  either  the  conductivity  or  the  thick¬ 
ness  of  the  crust  or  upper  lithosphere  where  the  electrical  conductivity  may 
be  lower.  This  is  due  in  large  part  to  the  band  limited  nature  of  the  seafloor 
electromagnetic  spectrum  and  the  probable  presence  of  ocean-generated  source 
fields  of  variable  but  unknown  spatial  morphology.  In  any  case,  the  modal 
structure  of  the  electromagnetic  field  induced  by  external  sources  is  the 
poloidal  magnetic  type  in  the  terminology  of  Chave  (1982)  and  is  quite  insen¬ 
sitive  to  low  conductivity  material. 

The  existence  of  low  conductivity  regions  in  the  upper  lithosphere  is 
suggested  by  laboratory  data  on  various  minerals.  Such  a  region  would  exert  a 
profound  influence  on  the  interpretation  of  both  motional  induction  fields 
caused  by  large  scale  oceanic  flows  and  magnetotelluric  experiments  in  the 
ocean-continent  transition  region  (Cox  1980,  1981).  If  the  conductivity  of 
any  possible  lithospheric  channel  were  low  enough  and  its  horizontal  extent 
were  basin-wide,  then  potential  military  applications  for  secure,  electromag¬ 
netic  pulse-proof  communications  would  arise. 

The  background  magnetic  spectral  density  in  the  frequency  band  centered 
-23  2 

near  1  Hz  is  about  10  (V/m)  /Hz  at  the  floor  of  the  deep  ocean  (Cox  ££  al. 
1978).  This  is  at  least  one  million  times  lower  than  is  observed  on  land  due 
to  the  shielding  effect  of  the  ocean  and  the  absence  of  any  significant 
seafloor  noise  sources.  This  led  Cox  (1980)  to  propose  the  use  of  seafloor- 
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based  controlled  electromagnetic  sources  using  frequencies  within  two  decades 
of  1  Hz  to  probe  the  conductivity  of  the  oceanic  crust  and  lithosphere.  Two 
approaches  using  vertical  (Edwards  e£  a_l  1981)  and  horizontal  (Young  and  Cox 
1981,  Chave  and  Cox  1982)  electric  dipole  (VED  and  HED)  sources  have  been  pro¬ 
posed  or  applied  to  date. 

In  this  report  the  theory  of  Chave  and  Cox  (1982)  is  extended  to  the  case 
where  either  the  source  or  receiver  is  located  in  the  oceanic  crust,  as  in  a 
borehole.  Numerical  calculations  are  presented  which  indicate  that  very  lit¬ 
tle  is  gained  by  the  use  of  boreholes  unless  low  conductivity  (  <  0.00S  S/m) 
material  is  penetrated  for  a  significant  distance,  an  unlikely  occurrence 
unless  the  hole  is  deeper  than  1500-2000  m.  In  particular,  a  VED  source 
located  in  a  borehole  with  receivers  placed  on  the  seafloor  yields  generally 
weaker  signals  than  HED  sources  and  receivers  both  located  on  the  seafloor.  A 
seafloor-based  HED  source  with  borehole  receiver  shows  greater  promise  in 
terms  of  signal  strength,  however  engineering  limitations  for  this  arrangement 
are  more  severe.  V"  ~ 


2.  Hwgrx 


The  relevant  theory  is  an  extension  of  that  presented  by  Chave  and  Cox 
(1982).  The  Maxwell  equations  of  electromagnetism  may  be  written  in  the  form 

V'B  -  0  (1) 

V  z  B  *  uoE  +  uJ°  (2) 

VxE  =-atB  (3) 

where  B  is  the  magnetic  induction,  E  is  the  induced  electric  field,  J®  is  the 
specified  or  impressed  source  electric  current,  fi  is  the  magnetic  permeabil¬ 
ity,  usually  taken  as  the  free  space  value  and  o  is  the  electrical 
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conductivity.  If  the  electrical  structure  of  the  earth  is  assumed  to  be  one 
dimensional,  then  (1)  -  (3)  can  be  rewritten  in  terns  of  two  scalar  functions 


II  end  f*  (Chave  1982) 

B  *  V  x  (IK)  +9x9  x(l*z)  (4) 

e  *  9ja„n/»ia  -  T/oj  -an*  -  9  x  (atf)  (5) 

I  z  |  t  t 

where 

7°  *  Jjz  +  9hT  +  9x  (Yz)  (6) 

and  the  governing  differential  equations  are 

V&I  +  4«a  (a  11/ no)  -  |io3.n  *  -  pi®  +  (T/o)  (7) 

n  z  z  z  z  z 

-  n«dt4*  *  -  (lY  (8) 

The  source  current  terns  in  (6)  are  the  solutions  of  the  Poisson  equations 

vJt  -  V^’J0  (9) 

v£y  -  -  (9x  J°)z  (10) 


The  separation  of  (4)  -  (10)  is  far  more  than  a  nathenatical  convenience.  The 
scalar  terms  II  and  I*  represent  the  two  independent  nodes  of  induced  current 
flow  in  a  one-dimensional  earth.  The  4*  function  represents  the  transverse 
electric  (TE)  or  poloidal  magnetic  mode  where  electric  current  flows  parallel 
to  the  earth's  surface  everywhere,  with  coupling  to  deep  structure  occurring 
by  induction.  The  I]  function  represents  the  transverse  magnetic  (TM)  or 
toroidal  magnetic  mode,  where  the  induced  electric  current  flows  in  loops  con¬ 
taining  the  vertical  direction.  TE  modes  are  analogous  to  the  magnetotellu- 
ric  situation,  and  are  quite  insensitive  to  low  conductivity  regions  since 
they  can  couple  across  them  by  mutual  induction.  TM  modes  are  analogous  to  DC 
currents,  and  are  by  contrast  very  sensitive  to  low  conductivity  zones  since 
vertical  electric  currents  cannot  penetrate  them. 
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For  artificial  soorce  work  there  are  four  fundamental  source  types:  the 
vertical  and  horizontal  electric  dipoles  (VED  and  HED)  and  the  vertical  and 
horizontal  magnetic  dipole  (VMD  and  HMD).  The  VED  and  HMD  can  induce  only  TM 
modes,  the  VMD  can  induce  only  TE  modes,  while  the  HED  produces  a  combination 
of  the  two.  Chave  and  Cox  (1982)  showed  that  the  HED  produces  the  best  cou¬ 
pling  between  ocean  and  earth,  while  Edwards  e£  al .  (1981)  considered  the  DC 
limit  and  a  VED  source.  The  VMD  and  HMD  sources  are  not  practical  in  the 
ocean  for  engineering  reasons,  except  possibly  on  a  very  small  scale.  In 
succeeding  sections,  a  VED  source  located  in  the  ocean  crust  with  seafloor 
receivers  and  an  HED  source  on  the  seafloor  with  borehole  receivers  will  be 
considered  in  detail. 


2.1  VED  Borehole  Source 

The  ocean  is  assumed  to  occupy  the  region  0.£z<H  with  conductivity  oQ,  air 
is  the  region  z>H,  while  the  VED  source  stretches  from  z  =  -L  to  z  =  0  with 
crustal  conductivity  o^,  and  an  arbitrary  conductivity  structure  is  located 
below  z  =  HL.  The  electromagnetic  fields  must  satisfy  (7)  since  only  TM  modes 
need  be  considered,  yielding 

-iuppoci  *  (id 

where  the  right  hand  side  vanishes  outside  -Li z<0  and  eiwt  sinusoidal  time 
dependence  is  assumed.  The  appropriate  solution  of  (4)  -  (6)  and  (11)  is  most 
easily  obtained  by  rewriting  the  Green  functions  obtained  by  Chave  (1982)  in 
the  form 

-Pjlz-z'l  -2p.L  -p.(z+z’)  P,(z+z')  -2{J..L  P-lz-z’l 

kim~ - ^n: - (12) 

^l^uV  > 

where  z  and  z'  are  the  observation  and  source  heights. 
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=  Jk2  +  iu,V.  (13) 
in  the  region  of  conductivity  o.  with  horizontal  wavenumber  k,  and  the  spatial 
dependence  is  recovered  from  the  Bessel  transform 


fl  *  /dkJ0(kp)kft  (14) 

0 

where  Jq  is  the  Bessel  function  of  the  first  kind  of  order  0  and  p  is  the  hor¬ 
izontal  range.  R^  and  R^  are  the  reflection  coefficients  at  the  z  =»  0  and 
z  *  -L  interfaces  respectively,  given  by 

atanhpQH-1 

atanhPQH+l 


h 

— K  -  1 

_  fl _ 

“l  mh 
-*%.  +  1 
°1 

where  a  ■  ^i°o^^O°l  ao<*  *  »ode  response  function  described  by  Chave 

and  Cox  (1982)  and  Chave  (1982)  which  contains  all  of  the  necessary  informa¬ 
tion  about  the  conductivity  below  z  ■  L.  For  a  finite  VED  extending  from 
z  =  -L  to  z  *  0  the  Green  function  (12)  is  integrated  to  yield 


n  -  77  i  i + 
4"  ip: 
i  1 


P.z  -2P-L  -p.z  -p  (L  +  z)  -p.  (L 

(Rp  -  l)(e  -t-  BLe  e  *)♦«--  l)  (c  +  Rne 


Pl(1  ‘Vl* 


-2H  L 


where  I  is  the  source  current.  This  expression  holds  only  in  the  source 
region  (-L  1  z  1  0),  and  the  solution  in  the  ocean  is  of  the  form 


A  “*0*  V 

fi  -  Ae  +  Be 


where  A  and  B  are  obtained  by  matching  the  boundary  conditions  at  z  ■  0  (con¬ 
tinuity  of  fl  and  3zn/o). 


After  tome  tedious  algebra,  the  results  are 


where 


i.Io! _ r-A 

B2  2P0l  2  ~2hL 

1  j»i(1  -  Vl*  > 

l.il _ _ 

2pol  1  -2^iL 

1  l»i<2  -  Vl«  > 


_  S0  -P.L1 

r  *  V®"  '  1>U  *  V  ’  *  ®L  •  1HRD  *  1>e  j 

Pj  '2P1L  "PlLl 

A  s  ^l(ru  -  1)(1  -  V  >  +  (kl  -  1)(Ru  -  1)6  J 

By  symmetry,  only  three  nonzero  electromagnetic  field  components  are 


obtained  from  (4)  and  (5) 


L-’< 


E,  •  4^  f0ia l“»>k  »o  1°*  -  j 

I  "  3  -P0X  ^0*1 

Ez  =  4no  ^dkJ0<kp)k  LCe  +  j 

»Qlmf  2  ~P0Z  V 

-  “/dkJ^kpU  jCe  +  De 

^O1  M0I 

where  C  *  -“pA  and  D  =  ~-B. 


2.2.  Baiahalfl.  Measurement  a  Seafloor  ID  Source 


The  mathematical  development  for  this  case  is  also  very  similar  to  the 
seafloor  source-receiver  case  treated  in  detail  by  Chave  and  Cox  (1982).  The 
solutions  for  a  borehole  occupying  -L  i  z  <  0  are  given  by 

Pi*  to  ~Pi*  “2PiL 

MnP  *  TO  -P0Z'  *  +  V«  • 

f  -  -rrdjd kjjkpn™*  - -1 - ttt —  (23) 


JTM  'pOz'  e 


Plz  TM  _P1Z  "2plL 

+  V  e  e 


1  -  Ru  \  e 


(24) 


for  an  infinitesimal  source  at  0  <  z'  <  ®  with  a  strength  p.  The  ex  ension  to 


a  finite  source  is  detailed  in  Appendix  A  of  Chave  and  Cox  (1982) 


TM 

Tg  are  seafloor  transmission  coefficients  given  by 


~TE  . 
TU  and 


„te 


2f*r 


P0  +  h 


(25) 


*0 

2— 


_TM 


D 


!n,!i 


TE  TM 

while  Rjj  and  are  the  seafloor  reflection  coefficients 


■? 


»1  -  >0 


Pi  +  pf 


(26) 


Pi  _  £o 


R 


TM 


!i  +  !o 


TM 

and  R^  and  R^  are  the  z  *  -L  reflection  coefficients 


Pi  A  -  l 


PXA  +  l 


(27) 


Pi 


K  -  1 

fl _ 


Pi 


— K  +  1 


with  the  response  functions  A  and  K  given  in  Chave  and  Cox  (1982).  The  six 


electromagmtic  field  components  follow  from  (4)  and  (S).  Defining  the  four 
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terms 


P1Z  TE 
6  *  \  6 


TE  ~P1Z  _2plL 


-V 


*  -  ^'2hL 


B(k)  =  - 


_TE  "Pi2  -2^1L 
~  \  6  e 

„TTLTE  "2P1L 


-<V' 


-C\n 


,  TW  -“l2  -2f>lL 

15,  «  » 

„TJLTM  _2P1L 


-V 


-«S\B 


P1Z  TM  “Plz  _2plL 
>  ~*L*  e 

*  - 


-V 


yields 


•  OS 

E  -  a  cos  A  /dkI(kp)T™p  D(k)  -  VdkJ  (kp^^Afk) 

p  [0  v  i  Pq  x  u  Pq 

2 

*  TF^l  1*  tw 

E .  =  a  sin  A  /dkI(kp)Ty  T±A(k)  -  VdkJ.tkpJT^VDtk) 

®  10  **0  1  u  1 

09 

E  ■  -a  cos  A  /dkJ,  <kp)T™k2C(k) 
z  Q  i  u 

B  -  b  sin  A  /dkI(kp)TjE^iB(k)  -  1  dkJ  <kp)  T™C(k) 

p  Lo  po  p  A  u 


=  b  cos 


A  |/dkI<kp)T™C<k)  -  “/dkJ.  (kp)T?^~B(k) ] 

LO  U  P0  1  °  P0  J 


B  -  b  sin  A  /dkJ.  (kpJT^J-AU)  (37) 

z  0  1  U  P0 

where  a  *  p/4 nc^,  b  =  Pqp/4jj,  A  is  the  azimuthal  angle  measured  with  respect 
to  the  source  wire,  and  I(kp)  *  kJQ(kp)  -  J^kpJ/p. 


2.1.  Nuaerical  Methods 
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The  field  expressions  (20)-(22)  and  ( 32 ) — ( 37 )  were  evaluated  nsing  Gauss 
quadrature  between  the  asymptotic  zero  crossings  of  the  Bessel  fnctions  JQ  and 
with  Pade  summation  of  the  resulting  series  of  partial  integrations.  For 
details,  see  Chave  and  Cox  (1982). 

3 .  Mode  1  Studies 

In  this  section,  the  electromagnetic  response  of  a  series  of  models  of 
the  upper  lithosphere  will  be  examined.  Similar  examples  to  those  used  by 
Chave  and  Cox  (1982)  will  be  used  when  possible  to  facilitate  comparison  of 
borehole  methods  with  the  seafloor  approach.  Some  additional  models  contain¬ 
ing  a  deep,  low  conductivity  channel,  as  postulated  for  the  continental  litho¬ 
sphere  by  Simmons  et  al..  (1980),  will  also  be  explored. 

Figure  1  shows  the  radial  and  vertical  electric  and  magnetic  fields  for  a 
seafloor  HED  source  and  receiver  as  a  function  of  range  at  a  frequency  of  1 
Hz.  The  ocean  model  has  a  conductivity  of  3.2  S/m  and  the  underlying  lithos¬ 
pheric  halfspace  has  a  conductivity  of  0.05  S/mj  the  corresponding  e-folding 

or  skin  depths  are  270  m  in  the  ocean  and  2.25  km  in  rock.  Note  that  the 

4  5 

infinitesimal  source  has  a  strength  of  1  A-m,  a  factor  of  10  -  10  smaller 

than  can  be  used  in  practical  experiments.  Attenuation  of  the  electromagnetic 
fields  in  a  quasi-exponential  fashion  with  range  is  apparent,  while  the  phase 
decays  at  about  one  radian  per  skin  depth  in  the  lithosphere.  The  direct 
field  propagating  in  the  highly  conductive  ocean  decays  very  rapidly  with 
range,  and  the  observed  fields  at  ranges  over  a  few  bn  are  those  traveling  in 
the  lower  conductivity  rock.  The  horizontal  components  of  both  the  electric 
and  magnetic  fields  are  at  least  ten  times  larger  than  their  corresponding 
vertical  components  at  all  ranges. 


Figure  2  shows  the  corresponding  results  for  a  borehole  VED  source  and 
seafloor-based  receivers.  The  source  occupies  the  upper  500  a  of  the  lower 
halfspace.  As  for  the  seaf loor-to-seaf loor  case,  the  horizontal  electric 
field  exceeds  the  vertical  component  by  about  a  factor  of  ten.  The  horizontal 
electric  and  magnetic  field  amplitudes  in  Figures  1  and  2  are  comparable, 
while  the  vertical  electric  field  is  significantly  stronger  at  range,  reflect¬ 
ing  the  slower  decay  rate  of  the  VED  field  (21),  when  compared  to  the  HED 
results  in  Chave  and  Cox  (1982).  Note  that  the  assumed  length  of  the  source 
is  dictated  by  practical  considerations  and  that  the  field  components  (20)- 
(22)  are  dependent  on  this  length. 

Figure  3  shows  the  results  obtained  by  borehole  measurement  of  a 
seafloor-based  HED  source.  By  contrast  to  the  seaf loor-to-seaf loor  and  VED 
cases,  the  vertical  electric  field  is  larger  than  the  horizontal  component 
near  the  source,  and  the  vertical  magnetic  field  is  comparable  to  its  horizon¬ 
tal  counterpart.  The  borehole  vertical  electric  field  is  smaller  than  the 
horizontal  component  of  Figure  1  at  all  ranges,  and  the  horizontal  magnetic 
components  are  more  nearly  equal  in  Figures  1  and  3.  Note  that  the  recipro¬ 
city  theorem  of  electromagnetism  requires  that  the  seafloor  horizontal  elec¬ 
tric  fields  produced  by  a  borehole  VED  source  and  the  borehole  vertical  elec¬ 
tric  field  produced  by  a  seafloor  HED  source  be  equivalent;  this  is  seen  in 
Figures  2  and  3. 

The  differences  observed  between  Figures  1-3  can  be  understood  in  a  qual¬ 
itative  fashion  by  a  simple  application  of  the  law  of  refraction.  For  a  hor¬ 
izontal  source  in  a  highly  conducting  region  like  the  ocean,  the  electric 
fields  are  largely  horizontal,  while  the  refracted  lines  in  the  lower  conduc¬ 
tivity  rock  will  be  nearly  vertical.  Similar  arguments  show  that  a  borehole 


VED  will  produce  largely  horizontal  fields  in  the  ocean.  This  qualitative 
argument  cannot  be  used  to  assess  the  relative  utility  of  the  three  cases 
treated  since  the  actual  field  magnitudes  depend  critically  on  the  electrical 
conductivity  structure  below  the  source. 

Chave  and  Cox  (1982)  presented  a  model  for  the  electrical  conductivity  of 
the  ocean  crust  based  on  ophiolite  and  laboratory  data.  The  seafloor-to- 
seafloor  electromagnetic  response  for  one  possible  model*  consisting  of  a  1  km 
pillow  basalt  layer  of  conductivity  0.05  S/m,  a  2.5  km  gabbro  layer  of  conduc¬ 
tivity  0.003  S/m,  and  a  2  km  sheeted  dike  region  of  conductivity  0.001  S/m 
overlying  a  0.0005  S/m  halfspace,  is  reproduced  as  Figure  4.  The  correspond¬ 
ing  result  for  a  borehole  VED  in  the  top  500  m  of  the  pillow  layer  is  shown  in 
Figure  5.  The  horizontal  electric  field  amplitude  is  in  all  instances  lower 
for  the  borehole  VED  source  at  long  range  and  comparable  at  a  few  km  range. 
The  relative  changes  in  both  amplitude  and  phase  produced  by  the  crustal  con¬ 
ductivity  structure  is  smaller  for  the  VED  source,  suggesting  decreased 
resolving  power.  In  particular,  the  very  large  phase  changes  of  Figure  4  are 
not  seen  in  Figure  5. 

Figure  6  displays  results  corresponding  to  Figure  4  for  the  azimuthal 
magnetic  field  for  a  seafloor  HED  to  seafloor  receiver  and  the  same  crustal 
model.  Very  similar  behavior  to  the  electric  field  result  of  Figure  4  is 
seen.  Figure  7  shows  the  same  field  component  from  the  same  source  at  a  depth 
of  500  m  in  a  borehole  in  the  upper  layer  of  the  crust.  The  magnetic  field 
amplitudes  are  similar  for  the  two  cases,  but  the  relative  changes  produced  by 
the  crustal  structure  are  very  much  smaller  for  the  borehole  receiver.  The 
large  amplitude  change  at  low  frequencies  for  the  seaf loor-to-seaf loor  path 
are  caused  by  the  geophysically  uninteresting  direct  wave  in  the  water. 


Figure  8  shows  three  models  of  the  upper  lithosphere  containing  a  subcru- 
stal  low  conductivity  channel.  In  all  cases  the  ocean  is  a  halfspace  of  con¬ 
ductivity  3.2  S/m,  while  the  crustal  model  of  Figure  4  is  simplified  to  a  S.5 
km  layer  of  conductivity  0.003  S/m.  This  neglects  the  high  conductivity  pil¬ 
low  basalt  layer  and  will  produce  an  optimistic  estimate  of  the  field  ampli¬ 
tudes.  Models  1  and  2  have  a  4  km  region  of  conductivity  0.0002  S/m  underly¬ 
ing  the  crust,  and  a  10  km  thick  low  conductivity  channel  at  10  **  S/m  for 
Model  1  and  10  ^  S/m  for  Model  2.  Model  3  differs  from  Model  2  in  the  subcru— 
stal  layer,  which  has  a  much  lower  conductivity  of  5  z  10  *  S/m.  This  value 
is  possible  only  if  mantle  rocks  contain  essentially  no  free  water. 

Figures  9-11  compare  the  amplitude  of  the  horizontal  electric  field  pro¬ 
duced  by  the  seafloor  HED  to  seafloor  receiver  with  the  borehole  VED  source  to 
seafloor  receiver  and  seafloor  HED  source  to  borehole  receiver  for  these  earth 
models.  In  each  case,  an  upper  500  m  crustal  conductivity  of  0.05  S/m  or 
0.005  S/m  is  used  for  the  borehole  cases;  the  former  is  consistent  with  meas¬ 
urements  from  DSDP  holes  in  the  upper  1  km  of  the  crust  (von  Herzen  e_£ 
1982),  while  the  latter  is  consistent  with  much  deeper  crustal  structure.  For 
Models  1  and  2  the  VED  source  produces  weaker  fields  than  the  seafloor  HED, 
while  this  situation  reverses  for  the  lower  subcrustal  conductivity  case  of 
Model  3.  The  borehole  receiver  and  seafloor  HED  consistently  produces 
slightly  higher  field  amplitudes,  especially  for  the  low  conductivity  crust, 
but  the  differences  are  less  than  a  factor  of  ten.  The  phases  (not  shown)  are 
consistent  with  the  results  of  Figure  4-7,  with  smaller  changes  occurring  when 
the  borehole  is  used. 


The  preceding  models  indicate  only  the  expected  signal  levels  and  rela¬ 
tive  changes  caused  by  conductive  structure.  The  separate  question  of  resolu¬ 
tion.  or  the  ability  to  quantitatively  estimate  conductive  structure  from 
field  measurements,  was  not  addressed.  The  fact  that  both  borehole  techniques 
produce  markedly  lower  relative  changes  in  the  electromagnetic  fields  when 
compared  to  the  seaflooar  method  suggests  that,  for  a  given  signal-to-noise 
ratio,  the  seafloor  technique  will  yield  better  results. 

The  estimation  of  various  experimental  parameters,  notably  antenna  resis¬ 
tances  and  required  source  voltage  and  current,  are  discussed  by  Cox  e_t  a  1 . 
(1980)  for  the  seafloor  HED  technique.  Some  simple  calculations  will  allow 
their  extension  to  borehole  methods.  For  the  borehole  VED,  the  transmitting 
antenna  is  located  in  a  medium  of  resistivity  P,  has  an  electrode  length  2L, 
and  a  radius  a,  and  approximating  its  shape  by  a  prolate  spheroid  yields  a 
resistance  to  infinity  of 


This  can  be  approximated,  when  a/L<<l,  as 

R  ■  4*  l0«  L^j  ”9) 

For  practical  purposes,  L  may  be  of  order  25  m,  while  a  is  the  radius  of  the 
borehole  (10  -  20  cm).  For  the  upper  oceanic  crust,  the  rock  resistivity  is 
20  0  -  m,  increasing  to  200  2  -  m  deeper  down.  This  yields  an  antenna  resis¬ 
tance  of  0.4  to  4  Q,  a  value  which  must  be  doubled  to  include  both  electrodes, 
and  the  intrinsic  resistance  of  the  antenna  wire  (typically  12)  must  be  added 
to  yield  a  total  resistance  of  2  -  10  2.  To  get  a  source  moment  of  105A-m, 
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similar  to  the  seafloor  QED,  will  require  a  source  current  of  200  A  for  a  500 
m  long  antenna.  This  in  turn  requires  a  generator  capable  of  delivering 
P  *  I2!!  watts,  yielding  P  *  80  -  400  kw.  Such  large  capacities  are  rarely 
available,  and  a  reduced  antenna  current,  and  concomitant  reduction  in  source 
moment,  will  undoubtedly  result.  It  should  be  noted  that  this  is  a  pessimistic 
estimate  since  L  was  taken  only  as  the  length  of  the  bared  wire  end  of  the 
antenna,  and  its  effective  length  will  be  larger  due  to  the  presence  of  con¬ 
ductive  seawater  in  the  borehole.  This  is  not  likely  to  change  the  results  by 
more  than  a  factor  of  2-4,  and  it  would  require  a  factor  of  over  20  to  make 

the  deep  crustal  source  practical.  Since  it  is  only  a  deep  crustal  VED, 

located  in  material  of  conductivity  below  .05  S/m,  which  yields  any  advantages 
over  the  seafloor  HED,  this  suggests  that  the  borehole  VED  source  is  not 
favored.  The  seafloor  HED  to  borehole  receiver  may  be  somewhat  more  practical 
if  deep  penetration,  into  the  lower  crust,  is  obtained.  Clearly  reception  of 
the  horizontal  electric  field  is  impractical  in  a  confined  hole,  and  the  vert¬ 
ical  electric  field  strength  is  not  large  at  long  source  receiver  ranges  from 
Fig.  3.  The  use  of  magnetometers  may  yield  better  results,  especially  if 
SQUID  types  or  the  collector-type  fluzgate  can  be  adapted.  The  motion  prob¬ 
lem,  which  is  a  serious  limitation  on  the  seafloor,  may  be  smaller  in  a 
borehole.  This  requires  further  study.  In  any  case,  only  deep  (  >  2  km  ) 
boreholes  in  the  ocean  crust  will  yield  any  real  advantage  over  the  seafloor- 

to-seafloor  method,  and  such  holes  are  hideously  expensive  to  drill.  For 

example,  DSDP  Hole  504B,  which  penetrated  to  1400.  m  in  the  ocean  crust  on  the 
Costa  Rica  Rift,  required  three  drilling  legs  of  1  1/2  -  2  mos.  each  on  Glomar 
Challenger . 
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Figure  1 


The  electric  and  magnetic  fielda  produced  at  the  seafloor  by  a  seafloor- 
bated  BED  source  as  a  function  of  range  for  a  source  frequency  of  1  Hz.  The 
radial  (solid  line)  and  vertical  (dashed  line)  components  are  shown  for  the 
electric  field  (left  panels)  at  an  azimuth  of  0°  and  the  magnetic  field  (right 
panels)  at  an  azimuth  of  90°.  Azimuth  is  the  angle  measured  with  respect  to 
the  end  of  the  source.  The  ocean  conductivity  is  taken  as  3.2  S/m,  while  the 
rock  halfspace  has  a  conductivity  of  0.03  S/m.  The  source  dipole  is  infini- 
tesimal  in  extent  with  a  source  moment  of  1  A-m.  (Taken  from  Chave  and  Cox, 
1982,  Figure  4.) 
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Figure  2 


T!ie  azimuthal  magnetic  field  (top),  radial  electric  field  (bottom,  solid 
line),  and  vertical  electric  field  (bottom,  dashed  line)  as  a  function  of 
range  at  a  frequency  of  1  Hz  produced  by  a  borehole  VED  source  and  observed  on 
the  seafloor.  The  source  is  assumed  to  occupy  the  upper  500  m  of  the  rock 
with  a  source  moment  of  1  A-m.  The  conductivity  model  is  that  of  Figure  1. 
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Figure  3 


The  electric  end  magnetic  fields  produced  at  S00  m  depth  in  the  rock 
half space  of  Figure  1  by  a  seafloor  HED  source.  The  top  panel  shows  the 
radial  (solid)  and  vertical  (dashed)  magnetic  fields  at  an  azimuth  of  90°, 
while  the  bottom  panel  shows  the  radial  (solid)  and  vertical  (dashed)  electric 
fields.  The  remaining  parameters  are  as  for  Figure  1. 
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Figure  4 


The  radial  electric  field  at  an  azimuth  of  0°  for  a  mature  oceanic  crust 
model  consisting  of  a  1  km,  0.0 5  S/m  pillow  and  flow  basalt  layer,  a  2.5  km, 
0.003  S/m  gabbro  layer,  and  a  2  km,  0.001  S/m  sheeted  dike  region  overlying  a 
0.0003  S/m  mantle  halfspace.  The  left  panels  show  the  fields  as  a  function  of 
source  frequency  at  ranges  of  4  km  (solid  line),  10  km  (short  dashed  line), 
and  16  km  (long  dashed  line).  The  right  panels  show  the  fields  at  ranges  of  3 
km  (short  dashed  line),  4  km  (solid  line),  5  km  (long  dashed  line),  and  8  km 
(dash-dot  line).  The  field  amplitudes  are  for  a  1  A-m  source  strength.  (Taken 
from  Chave  and  Cox  (1982),  Figs.  10  and  11) 
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Figure  5 


The  radial  electric  field  produced  at  the  seafloor  by  a  borehole  VED 
occupying  the  upper  500  m  of  the  same  model  of  Figure  4.  The  left  panels  shov 
the  field  at  ranges  of  4  km  (solid  line).  10  km  (short  dashed  line),  and  16  km 
(long  dashed  line).  The  right  panels  show  the  field  at  ranges  of  3  km  (short 
dashed  line).  4  km  (solid  line).  5  km  (long  dashed  line),  and  7  km  (dash-dot 
line).  The  field  amplitude  is  normalized  to  a  source  moment  of  1  A-m,  as  for 
Figure  4. 
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Figure  6 

The  azimuthal  magnetic  field  at  an  azimuth  of  0°  produced  at  the  seafloor 
by  a  seafloor  HED  source  for  the  same  crustal  model  as  in  Figure  4.  The  left 
panels  show  the  ranges  4  km  (solid  line),  10  km  (short  dashed  line),  and  16  km 
(long  dashed  line).  The  right  panels  show  the  ranges  3  km  (short  dashed 
line),  4  km  (solid  line),  S  km  (long  dashed  line),  and  7  km  (dash-dot  line). 
The  field  amplitude  is  normalized  to  a  source  moment  of  1  A-m. 
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Figure  8 


Three  models,  numbered  1-3  £rom  top  to  bottom,  for  an  oceanic  lithosphere 
containing  a  low  conductivity  channel.  For  the  high  and  low  conductivity 
borehole  cases  (see  text),  the  upper  S00  m  of  the  lithosphere  is  replaced  by  a 
layer  of  conductivity  0.05  S/m  and  0.005  S/m  respectively. 


Figure  9 


The  redial  electric  field  at  an  azimuth  of  0°  for  Model  1.  The  ranges  30 
km  (top)  and  70  ka  (bottom)  are  shown.  The  solid  line  represents  the  seafloor 
BED  to  seafloor  case,  while  the  borehole  VED  results  in  high  (0.03  S/m)  and 
low  (0.003  S/m)  conductivity  crust  are  shown  as  long  and  short  dashed  lines, 
and  the  borehole  BED  result  in  high  and  low  conductivity  crust  are  shown  as 
dash-dot  and  dash-dot-dot  lines  respectively.  The  field  amplitude  is  normal¬ 
ized  to  a  source  moment  of  1  A-m. 
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Figure  11 


The  radial  electric  field  at  an  azimuth  of  0°  for  Model  3.  See  Figure  9 
caption  for  details. 
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